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In an effort to understand the glass transition, the dy- 
namics of a non-randomly frustrated spin model has been 
analyzed. The phenomenology of the spin model is similar to 
that of a supercooled liquid undergoing the glass transition. 
The slow dynamics can be associated with the presence of ex- 
tended string-like structures which demarcate regions of fast 
spin flips. An entropy-vanishing transition, with the string 
density as the order parameter, is related to the observed 
glass transition in the spin model. 



The glass transition in supercooled liquids is heralded 
by anomalously slow relaxations with a time scale diverg- 
ing as the liquid freezes into the glassy state Q| . In recent 
years, there have been careful experimental and theoreti- 
cal studies aimed at understanding the structural aspects 
of this transition. The presence and nature of dynami- 
cal heterogeneities near the glass transition has been a 
dominant underlying theme of both simulations |^,^ and 
experiments [Q,^. Simulations in Lennard- Jones liquids 
1^ have shown the existence of string-like dynamical het- 
erogeneities and similar structures have been observed 
directly in colloidal glasses [||. In the Adam-Gibbs sce- 
nario, the glass transition is related to a phase transi- 
tion accompanied by the vanishing of configurational en- 
tropy 1^,0. An explicit connection between (a) dynam- 
ical heterogeneities, (b) anomalous relaxations and (c) 
the Adam-Gibbs scenario would provide useful insight 
into the nature of the glass transition. In this work, 
we present our analysis of a simple model where there 
are naturally occurring dynamical heterogeneities in the 
form of strings and where there is an entropy vanish- 
ing transition involving these structures. Monte Garlo 
simulations of the model show that there is a glass-like 
transition with diverging time scales and an anomalously 
broad relaxation spectrum. Analysis of the simulation re- 
sults provides strong evidence that the entropy- vanishing 
transition underlies the observed dynamical behavior. 

a. Model One of the simplest non-randomly frus- 
trated spin models is the triangular-lattice Ising anti- 
ferromagnet (TIAFM). The TIAFM has an exponen- 
tially large number of ground states and has a zero- 
temperature critical point P-jlO|]. The model studied 
in this letter is the compressible TIAFM (CTIAFM) 
in which the coupling of the spins to the elastic strain 
fields removes the exponential degeneracy of the ground- 



state. We solve the CTIAFM exactly within the ground- 
state ensemble of the TIAFM and show that there is 
an entropy-vanishing transition which involves extended 
structures. We then present results of simulations which 
indicate that the entropy-vanishing transition leads to 
glassy dynamics. 

The Hamiltonian of the CTIAFM is: 

H = Jy SiSj — ejy Ca} SiSj 
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Here J, the strength of the anti-ferromagnetic cou- 
pling, is modulated by the presence of the second term 
which defines a coupling between the spins and the ho- 
mogeneous strain fields eQ,a = 1,2,3, along the three 
nearest-neighbor directions on the triangular lattice. The 
last term stabilizes the unstrained lattice. The total 
number of spins in the system is given by N. The 
ground-state of the CTIAFM is a three-fold degener- 
ate striped phase where up and down spins alternate 
between rows and there is a shear distortion charac- 
terized by ei = e and 62 — — — e if the rows are 
along the direction of ei [ pTp^ ]. Within the manifold 
of the TIAFM ground-states, the competition between 
energy gained from the lattice distortions and the ex- 
tensive entropy of the TIAFM ground-states leads to an 
entropy- vanishing transition. The order parameter asso- 
ciated with this entropy-vanishing transition is the den- 
sity of extended string-like structures which characterize 
the TIAFM ground states. 

The string picture of the TIAFM ground states derives 
from a well-known mapping of these states to dimer cov- 
erings [[l3yi4[ . In a ground state, there is one unsatisfied 
bond per triangular plaquctte and the dimers are the 
fillcd-in bonds of the dual honeycomb lattice that cross 
the unsatisfied bonds of the triangular lattice, as shown 
in Fig.(^. Superposing a dimer configuration on a "stan- 
dard" dimer configuration where all the dimers are ver- 
tical [|l3|,0 leads to a string configuration (c/Fig.(^). 
Assuming spin- flip dynamics for the moment, the only 
spins that can be changed while the system remains in 
the ground-state manifold, are the ones which have a co- 
ordination of 3-3 (3 satisfied and 3 unsatisfied bonds). 
These are the fast spins in the system, and, as shown in 
Figure (|l|), are located at isolated kinks on the strings. 
The strings, therefore, play the role of dynamical hetero- 
geneities in this lattice model. 

b. Exact results We can solve the CTIAFM exactly 
within the restricted spin ensemble of the ground states 
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of the TIAFM. All the states in this ensemble can be 
classified according to the string density p = Ng/L where 
L is the linear dimension of the sample and is the 
number of strings [ p^ . The number of spin states {fl{p)) 
belonging to a particular string-density sector p has been 
shown to be exponentially large ^{p) — exp{N"f{p)), 
with N = L X L being the total number of spins. As 
shown in Fig. (||), the entropy 7(p) has a peak a,t p — 2/3 
In the CTIAFM (Eq. 0), the strain Cq couples to 
{SiSj)a- This average counts the number of strings along 
the a direction, i.e. the number of strings obtained by 
taking an overlap with the standard dimer state with 
all the dimers perpendicular to the a direction. Under 
periodic boundary conditions, only one of the three string 
densities is independent. 

+J - I + I + " -J + I + ' - ■ 

- [ + I© I + I - ^+ 
e] -\+ C? I - 1+ I -t^ 

-[^Isi - +1®]- + © 

S ^ ^ r ^ ^ 
e) - le I ® J - |e I e 

- I ^ - fe I a| - 1 + I © 
my — r (* I pfti — r f 




the coupling constant is increased, this minimum stays 
pinned at 2/3 and a second minimum starts developing 
at p = 0. The p = state stays metastable until at 
;i/Ti = (3/4) * 7(2/3) ~ 0.24 there is a first-order tran- 
sition from the p = 2/3 state to the p = state. The 
p = 2/3 state loses its metastability at a larger coupling 
given by /x/T* = V37r/12. At T*, the entropy vanishes 
and the order parameter shows a discontinuous change 
from p = 2/3top = 0. This transition is reminiscent of 
the transitions observed in p-spin spin glasses 




FIG. 2. (a) Entropy as a function of the string density from 
the work of Dhar et al.( |Q) (b)The dimensionless free energy 
/(p) for [1, — 0.18. T* for this value of fj. is 0.397. Temperature 
is measured in units of l/ks. 



FIG. 1. String representation of a TIAFM ground state. 
Dark bonds of the dual honeycomb lattice (c/text) are the 
dimers which divide unsatisfied pairs of spins. Light bonds 
define the "standard" dimer configuration. The strings are 
made up of these two types and extend across the system. 
The fast spins, the ones with 3 — 3 coordination, have been 
encircled. 

The strain field appears in the Hamiltonian as a purely 
Gaussian variable and can be integrated out, and the 
CTIAFM energy per spin, in the restricted ensemble, can 
be written in terms of the one, independent string density 
p: 



E{p) = -0,/2)iil-2pf + 2(1- pf)) 



(2) 



Here /z = e^J^/E. The energy function E{p) distin- 
guishes different string sectors and is minimized by p = 0. 
The entropy of the ground states, j{p), on the other hand 
favors the p = 2/3 sector and the competition between 
energy and entropy leads to the possibility of a phase 
transition. In the thermodynamic limit, the partition 
function Z = '^pexp~^^^P\ is dominated by the string 
density which minimizes f{p) = f3E{p) — j{p), where /3 
is the inverse temperature. The exact free energy corre- 
sponds to this minimum value of f{p) and the only rel- 
evant coupling constant in the problem is /3/i. For small 
values of the coupling constant, the function f{p) shown 
in Fig (||(b)), has only one minimum at p = 2/3. As 



These exact results show that in the CTIAFM, there 
is an entropy vanishing transition which defines the limit 
of stability of the homogeneous, liquid-like p = 2/3 state. 
In the Adam-Gibbs scenario, such a transition underlies 
the glass transition. The exact results are valid for the 
CTIAFM acting within the ground-state ensemble of the 
TIAFM. Defects [|l^, which correspond to triangles with 
all three bonds unsatisfied, can take the system out of 
the ground-state manifold. For low defect densities, how- 
ever, it is possible that the entropy-vanishing transition 
survives in some form and leads to slow glassy dynamics. 
We have investigated this scenario by performing Monte 
Carlo simulations of the CTIAFM. 

c. Simulations The parameters of the model were 
chosen to be J = 1, e = 0.6 and E ~ 2. These val- 
ues yield a small value of /i (= 0.18) and ensures that 
at the transition, T* {— 0.397), the defect density is 
low. The average defect number-density was measured 
to be ~ 0.04% at T = 0.45. We used Monte Carlo sim- 
ulations to study the dynamics of the supercooled state 
following instantaneous quenches to temperatures below 
Ti — 0.75. Spin-exchange kinetics was extended to in- 
clude moves which attempted changes of the strain fields 
Cq,. Details of the simulation algorithm have been pub- 
lished earlier ||l^ . System sizes ranged from 48x48 up to 
120x120. Unless otherwise stated, the results presented 
in this letter were obtained from 96x96 systems. 

As the glass transition is approached, global quantities 
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such as the energy per spin should exhibit anomalously 
slow relaxation processes. Fig shows the energy au- 
tocorrelation function CE{t,to) ~< E{to)E{t + to) >. 
For quench temperatures between T = 0.6 and T=0.47, 
CeCijio) is independent of the time origin to and has a 
stretched exponential form; exp"*^*/"^^^ . The stretching- 
exponent P decreases from 0.45 to 0.35 over this temper- 
ature range and the time scale te increases rapidly as is 
evident from the top panel of Fig (||). Below T ~ 0.47, 
the energy autocorrelation function depends on the time 
origin to, indicating that the equilibration times have be- 
come longer than our observation times. To illustrate 
the dependence on the waiting time to, we have shown 
the autocorrelation function averaged over three differ- 
ent ranges of to aX T = 0.45. The system is seen to relax 
more slowly for longer waiting times tg. This behavior 
of the energy autocorrelation function is similar to that 
of supercooled liquids and the temperature T = 0.47 is 
analogous to the laboratory glass transition temperature 
at which the equilibration time becomes longer than the 
observation time. In the CTIAFM, the proximity of this 
transition to T* suggests that the glassy dynamics is re- 
lated to the entropy vanishing transition. 
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FIG. 3. Top frame shows C'E{t,to) for three different tem- 
peratures where C_B(i,io) does not depend on to. The solid 
lines are stretched exponential fits. Bottom frame shows the 
waiting-time dependence of the energy autocorrelation func- 
tion at T = 0.45. The curves were obtained by averaging 
CE{t,to) over three different ranges of to- From bottom to 
top, these ranges are < to < 25000, 18000 < to < 48000 
and 50000 < io < 80000. 

The microscopic picture of the entropy- vanishing tran- 
sition is one where the string density vanishes. An anal- 
ysis of the string-density relaxation can, therefore pro- 
vide some insight into the nature of the dynamics at 



the glass transition. We find that the string-density au- 
tocorrelation functions are well described by exponen- 
tial relaxations, in contrast to the energy autocorrela- 
tion functions. Fig. (^) shows the results of our simu- 
lations for the relaxation times and fits to a power-law 
and a Vogel-Fulcher form Both fits yield a time- 

scale divergence at a temperature T ~ T*. The te ob- 
tained from the stretched exponential fits to CE{t,to) 
has a temperature dependence which tracks that of the 
string relaxation time. This observation suggests that 
the slow, non-exponential relaxations are a consequence 
of the freezing of the string-density relaxation which, in 
turn, is related to the entropy vanishing transition. The 
static susceptibility associated with the string density 
changes only by a factor ~ 2 over the temperature range 
in which the time scales change by a factor ~ 40, indi- 
cating that the entropy-vanishing transition suffers from 
anomalously strong critical slowing down. 
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FIG. 4. Temperature dependence of the string-density re- 
laxation time. The simulation results are shown with error 
bars. The solid line is a fit to the Vogel-Fulcher form and the 
dashed line to a power-law. 

In order to investigate the nature of the string- 
relaxations further, we measured the distribution, 
P{Ap), of Ap ~ p{t + to) — p{to), the deviation of the 
string density in the time interval t (Fig. (^). The 
distributions were generated by choosing different time 
origins to- The most striking feature of the distributions, 
observed at temperatures close to T = T*, is its non- 
Gaussian nature at intermediate times. At T = 0.55, the 
non-Gaussian feature is most pronounced at t = 4000. 
Beyond this time the distribution relaxes towards a Gaus- 
sian and for t > 8000, the distribution is time indepen- 
dent. At T = 0.47, a time-independent behavior is not 
observed for times as long as 30,000 and the distribu- 
tions are non-Gaussian at all intermediate times. In usual 
critical-point dynamics |p^ , one would expect to find a 
distribution of Ap which takes longer to reach its time- 
independent form as the critical point is approached and 
to find non-Gaussian behavior (within the limits of finite- 
size cutoffs) in the stationary distribution. In contrast, 
we observe the most pronounced non-Gaussian features 
at intermediate times. Drawing an analogy with critical 
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phenomena, this observation leads us to speculate that 
there exists a time-dependent length scale, ^(r) has a 
peak at a time tq{T). As T — > T*, the time scale tq and 
the height of the peak, £,{tq) appear to diverge. The exact 
nature of the divergence is difficult to extract from the 
current data. These apparent divergences indicate that 
the thermodynamic transition present in the zero-defect 
sector has been replaced by a dynamical transition [|l9|. 
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FIG. 5. Distribution of string-density deviation Ap for dif- 
ferent time intervals, t, shown at T = 0.55 and T — 0.47. The 
areas under the curves have been normalized to unity. 

d. Connection to real glasses In conclusion, our study 
of the CTIAFM provides strong indication that the slow, 
glassy dynamics in this model is associated with an 
entropy-vanishing transition involving extended, string- 
like structures. These structures are a manifestation of 
the frustration embodied in the nearest-neighbor, anti- 
ferromagnetic interactions, and the entropy-vanishing 
transition is a consequence of coupling to another degree 
of freedom, the lattice strain, which tends to remove the 
frustration in the system. The strings are naturally oc- 
curring dynamical heterogeneities since they demarcate 
regions of fast spin flips. If the relation between frustra- 
tion and dynamical heterogeneities is a generic feature of 
glass formers, then our observations would suggest that 
the clusters of mobile particles observed in Lennard-Jones 
simulations |^ and in colloidal systems ||^ should con- 
sist of the most frustrated particles in the system. In 
a Lennard-Jones mixture, these should be the particles 
with the least number of unlike bonds (if unlike bonds are 
energetically preferred) and in a colloidal system, these 
should be the particles which have coordinations that are 
furthest from being icosahedral. An experimental verifi- 
cation of this correlation between geometry and mobil- 
ity would be a direct test of the connection between the 
glass-transition and an entropy-vanishing transition in- 



volving extended structures forced in by frustration. 
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